0254-6124/2015/35(3)-350-06 Chin. J. Space Sci 空间 科学 学 报 


Wang Dazhao, Bai Weihua, Sun Yuediang, Meng Xiangguang, Wang Dongwei, Li Wei. Inversion algorithms and precision analysis 
of ROPP (in Chinese). Chin. J. Space Sci., 2015, 35(3):350-355, doi:10.11728/cjss2015.03.350 


ROPP 反 演 软件 算法 及 其 精度 分 析 
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摘 要 介绍 了 ROPP 反 演 软件 中 无 线 电 掩 星 反 演 的 算法 与 精度 分 析 . 采用 COSMIC 卫星 2008 年 1 月 1 日 
全 天 的 附加 相位 数据 , 反 演 得 到 折射 率 、 温度 、 压 强 与 湿度 等 参数 , 并 与 CDAAC 相应 结果 进行 对 比 . 实验 结果 表 
明 , 在 30 km 高 度 以 下 , 折射 率 、 压 强 和 湿度 的 相对 误差 在 2% 以 内 , 温度 误差 不 超过 2 K. 

关键 词 无线 电 掩 星 , 折射 率 , 大 气 参 数 , ROPP 

中 图 分 类 号 P128.14, V 448 


Inversion Algorithms and Precision 
Analysis of ROPP 


WANG Dazhaol” BAI Weihua2 SUN Yueqiang” 
MENG Xiangguang” WANG Dongwei2 LI Wei? 


1(University of Chinese Academy of Sciences, Beijing 100049) 


2(Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing 100190) 


Abstract The unique combination of global coverage, high precision, high vertical resolution, long- 
term stability and all-weather viewing of radio occultation will enhance the data sets of numerical 
weather prediction and improve the precision of Numeric Weather Prediction (NWP). And data 
assimilation requires preprocessing of the raw radio signals into the bending angles and refractivity. 
Then the atmospheric parameters such as temperature, pressure and humidity can be retrieved 
through data variational assimilation. The procedures for processing of radio occultation data, 
especially several important steps are introduced. And comparative analysis between these results 
and those from CDAAC shows that below the altitude of 30 km, the relative errors of refractivity, 
pressure and humidity are within 2% and the deviation of temperature is below 2 K. The results are 
also compared with the data from ERA-interim model of ECMWF and similar conclusion is obtained. 
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0 引言 


美国 喷气 推进 实验 室 (Jet Propulsion Laborato- 
ry, JPL) 在 19 世纪 60 年 代 首 次 提出 将 无 线 电 掩 星 
技术 应 用 于 地 球 大 气 主动 探测 趾 , 由 于 当时 技术 和 
经 济 条 件 的 限制 , 计划 暂时 搁浅 . 在 美国 与 前 苏联 分 
别 部 署 了 各 自 的 全 球 定位 系统 GPS 和 GLONASS 之 
后 , 利用 手 星 技术 探测 地 球 大 气 剖 面 成 为 可 能 .1993 
年 美国 UCAR (University Corporation for Atmo- 
spheric Research), UA (University of Arizonay) 与 JPL 
联合 制定 了 GPS/MET 计划 . 其 第 一 颗 实 验 卫星 Mi- 
croLabl 于 1995 年 成 功 发 射 . 通过 两 年 的 运行 , 证 明 
了 无 线 电 掩 星 技术 可 在 地 球 大 气 探测 与 研究 中 起 到 
重要 作用 7- 引 . 由 于 GPS/MET 计划 的 成 功 , NASA 
在 1995 年 资助 了 两 个 国际 计划 , 分 别 是 丹麦 Orsted 
卫星 和 南非 Sunsat 卫星 内 . 此 后 国际 上 较 出 名 的 无 
线 电 捧 星 卫 星 有 SAC-C, CHAMP, GRACE, Metop 
系列 和 COSMIC 等 . 值得 一 提 的 是 , 2013 年 中 国 成 
功 发 射 了 风云 三 号 卫星 , 首次 兼容 GPS 和 中 国 北斗 
导航 卫星 系统 开展 无 线 电 掩 星 探 测 . 

无 线 电 捧 星 具 有 全 球 覆 盖 、 高 精度 、 高 垂直 分 辩 
率 、 长 期 稳定 性 及 全 天 候 观测 等 优势 , 每 颗 LEO 卫 
星 每 天 大 约 可 以 获取 500 条 剖面 , 这 为 气候 研究 与 
数值 天 气 预 报 提供 了 大 量 数据 源 ， 目 前 无 线 电 掩 星 
研究 的 热点 之 一 是 , 如 何 将 掩 星 观测 到 的 数据 应 用 于 
现行 数值 天 气 预 报 系统 ， 其 关键 是 无 线 电 掩 星 反 演 
的 精度 应 满足 数值 天 气 预 报 系统 的 精度 要 求 . 目前 
比较 有 影响 力 的 GNSS 气象 学 研究 机 构 主 要 有 GFZ 
的 ISDC 数据 系统 中 心 、JPL 的 GENESIS 数据 系统 
中 心 及 UCAR 的 CDAAC 数据 系统 中 心 等 . UCAR 
与 中 国 台 湾 联合 发 射 的 COSMIC 任务 是 目前 公认 的 
在 无 线 电 掩 星 探 测 领域 取得 较 好 突破 的 任务 . 

ROPP 由 EUMESAT 下 属 科 研 机 构 GRAS 开 
发 , 其 最 初 目的 是 用 于 处 理 Metop 卫星 数据 , 但 在 实 
际 使 用 中 , 通过 调整 软件 配置 , 可 用 于 处 理 其 他 GPS- 
LEO 掩 星 数 据 ， 在 掩 星 数据 处 理 过 程 中 , ROPP 可 
以 对 无 线 电 掩 星 数 据 进行 预 处 理 和 质量 控制 , 还 可 
将 掩 星 数 据 同 化 到 NWP (Numeric Weather Predic- 
tion) 或 者 其 他 气候 模式 . ROPP 是 基于 Fortran 90 
的 开源 代码 , 在 Linux 系统 下 运行 是. 文献 [6] 介绍 
了 ROPP 各 模块 之 间 的 关系 , 以 及 其 中 主要 模块 的 
算法 原理 , 使 用 ROPP 处 理 了 一 组 COSMIC 数据 , 并 
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对 其 结果 进行 了 精度 分 析 . 在 文献 [6] 基础 上 , 根据 
实际 需求 , 本 文 将 ROPP 处 理 软件 移植 到 Windows 
中 并 用 Matlab 进行 编译 . 通过 以 上 调整 , 可 以 实现 
在 Windows 下 对 一 天 的 无 线 电 掩 星 数 据 进 行 批量 处 
理 , 共计 2354 次 掩 星 上 事件. 

使 用 ROPP 的 预 处 理 模 块 对 2009 年 1 月 1 日 
的 COSMIC 数据 进行 处 理 . 将 ROPP 的 反 演 结果 
与 OCC 反 演 结果 对 比 , 可 以 看 到 最 优 弯曲 角 、 折 射 
率 和 干 温 在 45 km 高 度 以 下 保持 了 良好 的 一 致 性 , 平 
均 相 对 偏差 在 0.4% 以 内 ; 而 在 45 km 高 度 以 上 误差 
较 大 , 其 原因 在 于 在 电离 层 修正 过 程 中 选择 了 不 同 的 
最 优 拟 合 MSIS 廓 线 . 同时 对 于 同一 天 的 数据 , 还 对 
同一 位 置 的 ECMWEF 分 析 数 据 进 行 了 处 理 , 得 到 折 

通过 对 比 ROPP 反 演 的 折射 率 和 ECMWEF 的 折 
射 率 , 可 以 看 出 在 多 路 径 效 应 发 生 区 域 , 相对 误差 可 
达 2% 以 上 ,和 干 温 误 差 更 为 明显 , 但 是 在 10~20km 高 
度 内 , 误差 较 小 , 不 超过 1K. 本 文 使 用 的 数据 为 2008 
年 1 月 1 日 的 COSMIC 数据 , 采用 的 ROPP 版 本 
为 v6.1, 主要 进行 了 ROPP 反 演 结果 与 CDAAC 数 
据 的 对 比 , 以 检验 ROPP 与 CDAAC 数据 处 理 流 程 
的 异同 及 其 精度 . 


1 无 线 电 掩 星 数 据 处 理 流程 


从 GNSS 发 射出 的 信号 在 穿越 电离 层 和 大 气 层 
后 , 被 LEO 接收 机 接收 . 由 于 电离 层 与 大 气 层 的 
折射 与 衍射 效应 导致 GNSS 信和 号 出 现 相位 延迟 . 
在 25km 高 度 以 上 , 传统 几何 光学 反 演算 法 根据 观 
测 到 的 相位 延迟 和 卫星 (GNSS 与 LEO) 相对 位 置 
与 速度 , 可 以 得 到 信和 号 的 弯曲 角 . 但 是 在 25 km 高 度 
以 下 , 由 于 水 气 含量 较 高 , 出 现 了 多 路 径 效 应 和 超 折 
射 , 几何 光学 反 演 不 适用 , 这 时 需要 进行 波动 光学 反 
演 , 从 而 得 到 25 km 高 度 以 下 的 弯曲 角 痢 面 . 由 于 电 
离 层 的 影响 , 在 反 演 大 气 参数 时 要 进行 电离 层 修正 . 
将 Ll 和 L2 波段 的 弯曲 角 插值 到 同一 碰撞 参数 层 之 
后 , 进行 线性 组 合 与 统计 最 优化 , 可 得 到 最 优 弯 曲 角 ， 
即 消 除了 电离 层 影响 的 大 气 弯曲 角 . 利用 Abel 逆 变 
换 , 可 以 得 到 折射 率 剖 面 . 由 于 折射 率 受到 温度 与 压 
强 的 影响 , 在 没有 水 气 的 情况 下 , 即 在 15km 以 上 高 
度 时 , 折射 率 仅 受 到 干 温 和 干 压 的 影响 , 结合 流体 更 
力学 方程 , 可 以 得 到 15km 高 度 以 上 的 干 温 和 干 压 . 
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但 是 在 15km 高 度 以 下 , 水 气 的 影响 是 不 能 忽视 的 ， 
仅 由 两 个 方程 无 法 得 到 温度 、 干 压 和 湿 压 三 个 参数 ， 
这 时 需要 使 用 大 气 模 式 的 背景 数据 , 对 掩 星 反 演 结果 
进行 变 分 同化 . 通过 积分 算 子 , 实现 同化 表示 成 分 析 
变量 函数 的 任意 观测 量 .，ROPP 中 可 以 对 弯曲 角 和 
折射 率 进行 同化 , 在 此 选择 对 大 气 折射 率 进行 同化 ， 
因为 对 折射 率 进行 同化 所 需 附 加 资料 仅 为 点 值 坐标 ， 
观测 算 子 相对 简单 , 容易 实现 . 

ROPP 包括 输入 /输出 模块 fopp-io)、 预 处 理 模 
块 (ropp-pp)、 一 维和 二 维 前 向 模块 (ropp-fm)、 质 
量 控制 工具 (ropp-qc) 和 一 维 变 分 反 演 模块 (rop- 
p-ldvar). 在 附加 模块 即 工具 模块 (ropp-utils) 中 
有 ROPP 模块 的 部 分 或 全 部 实用 程序 . 图 1 给 出 
了 ROPP 各 模块 间 的 相互 关系 以 及 各 模块 的 主要 画 
数 . 其 中 主要 的 是 ropp-pp 模块 和 ropp-ldvar 模块 . 
ropp-pp 模块 的 主要 输入 为 LEO 和 GPS 卫星 的 位 
置 和 速度 及 附加 相位 , 输出 主要 是 折射 率 、 干 大 气温 
度 和 压强 .ropp-ldvar 模块 的 主要 输入 包括 观测 数 
据 和 背景 数据 , 输出 为 湿 大 气温 度 、 压 强 和 湿度 等 . 
1.1 ”数据 处 理 流程 

无 线 电 掩 星 数 据 处 理 流程 如 图 2 所 示 . 可 以 看 
出 , 由 LEO 接收 机 接收 到 的 附加 相位 数据 , 反 演 得 到 
弯曲 角 和 折射 率 剖面 , 然后 得 到 没有 水 气 分 量 的 干 温 
和 干 压 . 在 得 到 折射 率 剖 面 后 , 对 折射 率 剖 面 进行 一 
维 变 分 , 进而 得 到 带 有 水 气 分 量 的 温度 、 压 强 及 湿度 


信息 . 


ropp_pp 
ropp_pp_occ tool 
ropp_pp_invert tool 
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1.2 主要 步 又 算法 说 明 

(1) 特定 任务 修正 

CHAMP L2 数据 有 很 高 的 仪 带 噪声 , 因此 L2 振 
幅 用 Ll 代替 . 光线 高 度 在 8~11km 以 下 时 , COS- 
MIC 和 GRAS 采用 了 开 环 技术 上 且 分 别 有 自 己 的 原始 
采样 模式 ， 该 区 域内 的 附加 相位 测量 反 演 要 去 除 在 
测量 过 程 中 加 到 每 个 数据 样本 中 的 伪 随 机 数 0 或 x. 
这 个 过 程 可 以 使 用 内 部 衍生 的 导航 数据 或 参考 导航 
数据 的 外 部 来 源 . 

(2) 数据 截断 

数据 截断 用 来 剔除 信号 跟踪 误差 非常 大 的 数据 . 
截止 的 标准 基于 测量 相位 和 通过 平滑 的 几何 光学 计 
算出 来 的 弯曲 角 和 碰撞 参数 剖面 的 大 概 佑 计 . 别 除 从 
振幅 低 于 最 大 振幅 的 某 一 比例 (例如 0.1) 时 刻 开 始 
的 数据 .剔除 从 平滑 弯曲 角 剖 面 大 于 特定 碰 拉 高度 
和 直线 切线 高 度 达 到 临界 值 (例如 0.05) 时 刻 开始 的 

(3) 无 线 电 全 息 滤波 

在 进行 几何 光学 反 演 前 , 应 对 无 线 电 掩 星 数 据 进 
行 处 理 , 将 数据 中 的 噪声 降低 并 去 除 较 差 的 数据 样 
本 . 在 此 ROPP 使 用 了 一 种 基于 正则 变换 (CT2) 和 
无 线 电 全 息 聚 焦 合成 孔径 法 的 无 线 电 全 息 法 滤波 算 
法 . 首先 将 接收 到 的 信号 进行 CT2 处 理 , 将 转换 空 
间 中 的 信号 除 以 参考 信号 , 然后 与 一 个 傅 里 叶 滤波 窗 
进行 卷 积 , 再 乘 以 参考 信号 从 而 重 构 初始 信号 . 该 算 
法 很 容易 实施 , 而 且 在 数值 模拟 中 得 到 了 验证 门 . 


preprocess module 


IO module 


utils module 


ld-VAR module 


ropp_fm_bg2ro_1d 
ropp_fm bg2ro 2d 


ropp_ldvar % 


ropp_ldvar refrac 


ropp_ldvar_bangle 


| | 
ropp2ropp 
ropp2bufr| bufr2ropp 


forward module 


图 1 ROPP 模块 及 其 相互 关系 
Fig.1 Relationships of the modules of ROPP 
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data input 


excess phase, 
SNR, 
position and velocity 
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compute the occultation 
point 


Specific correction 


statistical optimization 


ionospheric correction and optimization 
bending angles 


data cutoff 


Wl 


Abel inverse integral 一 refraction ; 


radio holographic filtering 
of L1 and L2 


hydrostatic integral 


dry temperature, 
and dry pressure, 


图 2 无线 电 掩 星 数据 处 理 流程 


Fig.2 Data processing procedures of the radio occultation 


(4) 几何 光学 

在 水 气 含量 较 少 的 海拔 高 度 上 , 没有 出 现 多 路 径 
效应 时 , 弯曲 角 是 碰撞 参数 的 单 值 函数 . 在 折射 率 的 
局 部 球 对 称 假设 下 , 结合 GNSS 和 LEO 卫星 的 相对 
位 置 和 速度 , 根据 Snell 方程 , 最 终 可 以 得 到 弯曲 角 
与 碰撞 参数 加， 传统 的 几何 光学 反 演 在 水 气 含量 较 
低 的 高 度 范围 内 , 反 演 精度 较 高 且 算 法 实现 容易 . 

(5) 波动 光学 

当 水 气 含量 变 高 时 , 产生 多 路 径 效 应 和 超 折射. 
此 时 单 信号 假设 及 瞬时 频率 与 观测 的 对 应 性 均 不 成 
立 , 即 观测 弯曲 角 与 碰撞 参数 不 是 一 一 对 应 的 . 基于 
波动 方程 的 波动 光学 反 演 方法 就 是 为 克服 上 述 多 路 
径 效 应 而 设计 的 . 主要 的 波动 光学 反 演算 法 有 衍射 
法 名、 后 传播 反 演 、 滑 动 频谱 法 、 全 谱 反 演 和 正则 
变换 回 . 这 五 种 算法 中 , 全 谱 反 演 与 正则 变换 的 反 演 
结果 最 好 , 而 全 谱 反 演 的 计算 效率 略 高 于 正则 变换 . 
在 此 选择 了 第 二 类 正则 变换 算法 . 

(6) 电离 层 修正 与 统计 最 优化 

电离 层 修正 残余 误差 与 接收 机 噪声 是 无 线 电气 
星 的 一 个 主要 噪声 源 , 这 些 噪声 会 影响 35~40km 高 
度 范 围 内 的 反 演 精度 . ROPP 使 用 的 电离 层 修 正 与 减 
小 噪声 的 算法 包括 以 下 细节 : 对 气候 模式 进行 拟 合 ; 


使 用 50km 以 上 的 无 线 电 掩 星 信 号 对 信和 号 与 噪声 的 
协 方差 矩阵 进行 动态 估计 ; 根据 背景 弯曲 角 、 信 和 号 与 
噪声 的 协 方差 矩阵 , 对 中 性 大 气 和 电离 层 弯 曲 角 进行 
统计 最 优化 分 析 09. 


2 无 线 电 掩 星 反 演 实例 


本 次 实验 选取 了 COSMIC 卫星 2008 年 1 月 1 
日 全 部 共 2354 次 大 气 掩 星 事 件 ， 其 中 下 降 掩 星 次 
数 为 967， 上 升 掩 星 次 数 为 1387. 探测 全 程 高 度 
为 0~60km, 对 COSMIC 数据 采用 了 开 环 跟踪 处 理 ， 
CDAAC 数据 处 理 成 功率 为 75.8%. 

从 包含 附加 相位 信息 的 文件 中 反 演 出 干 温 、 干 
压 . 再 使 用 背景 数据 进行 一 维 变 分 , 得 到 湿 温 、 湿 压 
与 湿度 信息 .将 反 演 结果 与 经 过 CDAAC 数据 处 理 
得 到 的 温 湿 压 进行 对 比 , 即 可 得 到 反 演 结果 的 平均 偏 
差 和 标准 差 , 结果 如 图 3 所 示 . 

图 3 中 红色 曲线 表示 ROPP 反 演 结果 与 C- 
DAAC 的 平均 偏差 , 其 左右 两 条 蓝 色 曲 线 分 别 表示 
平均 偏差 减 去 标准 差 、 平 均 偏 差 加 上 标准 差 的 值 . 从 
图 3 可 以 看 出 , 在 0~45km 高 度 折 射 率 反 演 结果 较 
好 , 折射 率 相对 误差 的 平均 偏差 保持 在 1% 以 内 . 
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Fig.3 Comparison of dry atmosphere inversion 


而 干 温 误 差 的 平均 偏差 在 0~35 km 高 度 范围 内 保持 
在 1K 之 内 , 干 压 相对 误差 则 保持 在 1% 之 内 . 

从 文献 [6] 可 以 看 出 , 对 于 单条 廓 线 , ROPP 反 
演 结果 与 CDAAC 的 结果 对 比 显示 , 在 40 km 高 度 以 
下 , 折射 率 相 对 误差 可 以 保持 在 1% 之 内 , 而 在 40 km 
以 上 , 其 迅速 增 大 , 可 以 达到 一 3%. 
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在 数据 处 理 过 程 中 , 主要 出 现 异常 的 问题 在 于 ， 
在 对 ROPP 反 演 的 数据 进行 一 维 变 分 同化 时 , 由 于 
质量 控制 参数 的 不 一 致 , 导致 ROPP 与 CDAAC 对 
于 同一 条 亡 线 的 质量 产生 了 不 同 的 结果 , 在 观测 数据 
与 背景 数据 进行 匹配 时 出 现 异 常 . 

随后 使 用 ECMWEF 业务 化 预报 的 模式 的 数据 作 
为 背景 数据 , 对 折射 率 进 行 一 维 变 分 同化 , 得 到 了 包 
含有 湿度 信息 的 温度 、 压强 和 湿度 , 结果 如 图 4 所 示 . 

由 图 4 可 以 看 出 , 湿 大 和 气 的 反 演 结果 要 略 差 于 
干 大 气 反 演 结 果 . 但 是 在 35 km 海拔 高 度 内 , 温度 误 
差 可 以 控制 在 2K 内 , 而 压强 与 湿度 的 相对 误差 均 小 
于 2%. 此 次 实验 中 还 将 反 演 结果 与 ECMWEF 再 分 
析 模 式 ERA-interim 的 数据 进行 了 对 比 , 得 到 的 结果 
类 似 . 在 文献 [6] 中 , 对 于 多 路 径 效应 发 生 的 区 域 , 温 
度 误差 较 大 , 可 达到 2K, 而 此 次 实验 中 , 对 于 多 路 径 
效应 发 生 区 域 , 误差 得 到 了 明显 改善 . 


3 结语 


通过 反 演 统计 分 析 可 知 , 折射 率 廊 线 和 干 大 气温 
度 、 压 强 的 反 演 精度 在 0~35 km 内 较 高 , 折射 率 相 
对 误差 在 1% 之 内 , 压强 相对 误差 在 2% 之 内 , 干 温 
误差 在 1K 左右 . 而 在 35km 高 度 以 上 , 反 演 精度 
有 所 降低 . 湿 大 气 反 演 结果 表明 , 在 10km 以 下 , 温 
度 误差 在 0.5~1.5K 以 内 ; 在 10~25 km 高 度 , 温度 
误差 不 超过 1K; 在 25 km 高 度 以 上 , 温度 误差 增 大 ; 
在 35km 高 度 以 下 , 压强 相对 误差 在 2% 之 内 . 比 湿 
误差 也 在 2% 之 内 . 

无 线 电 掩 星 反 演 在 30km 高 度 以 上 的 主要 
误差 来 源 是 电离 层 修正 不 完全 和 接收 机 跟踪 误差 ， 
而 10km 高 度 以 下 的 误差 主要 是 由 于 大 气 水 气 模糊 
度 无 法 确定 . 而 大 气 参 数 (温度 、 压 强 和 比 湿 ) 误差 
比 折射 率 误 差 大 , 其 主要 原因 是 在 数据 同化 过 程 中 采 
用 了 背景 数据 , 引入 背景 数据 误差 . 背景 数据 误差 对 
比 湿 的 影响 较 小 , 而 对 温度 和 压强 的 影响 较 大 . 比 湿 
在 10km 高 度 以 上 几乎 为 0, 而 在 10km 高 度 以 下 因 
水 气 梯度 变化 较 大 , 导致 误差 发 生 了 剧烈 波动 . 主要 
误差 来 源 包括 电离 层 修正 的 影响 、 接 收 机 噪声 、 水 气 

由 上 述 实验 可 以 看 出 , ROPP 反 演 软件 的 预 处 理 
模块 对 无 线 电 掩 星 数 据 的 处 理 结果 与 CDAAC 提供 
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Fig.4 Comparison of wet atmosphere inversion 


的 结果 较为 接近 , 而 一 维 别 . 在 一 维 变 分 过 程 中 出 现 
误差 的 可 能 原因 是 , 所 使 变 分 模块 结果 有 了 一 定 差 用 
的 背景 数据 不 一 致 或 者 同化 算 子 的 选择 不 一 致 . 在 此 
选择 的 背景 数据 分 层 仅 为 60 层 , 而 CDAAC 反 演 出 
来 的 数据 分 层 为 400 层 , 其 垂直 分 辨 紊 明显 优 于 本 文 
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所 选择 的 背景 数据 . 另外 , 选择 弯曲 角 同 化 与 折射 率 
同化 , 两 者 差别 也 会 造成 反 演 结果 的 不 同 . 

在 对 ROPP 反 演 软件 进行 系统 移植 与 部 分 修改 
之 后 , 还 可 以 根据 掩 星 事件 发 生 的 经 纬度 和 时 间 分 别 
进行 统计 分 析 . 可 以 量化 描述 赤道 地 区 、 温带 与 极地 
的 反 演 结果 误差 并 进行 对 比 , 也 可 量化 描述 同一 地 区 
不 同 地 方 时 的 误差 . 进一步 将 利用 一 个 月 左右 的 数 
据 , 对 不 同 纬度 带 的 反 演 结果 进行 分 析 . 
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